UNCLASSIFIED 


AD  401188 


if  iJm 

DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIEN  IFIC  AND  TECHNICAL  INFORMATION 

CA  MERON  STATION.  ALEXANDRIA.  VIRGINIA 


UNCLASSIFIED 


mnci:  Vban  gorermmt  or  othar  drawiagB,  ■paei- 
fleatiaoa  or  ottaar  data  ara  uaad  for  any  piupoaa 
othar  than  in  eonnaetian  with  a  daflnltaly  ralatad 
goranaant  proeunaant  oparation,  tha  U.  8. 
Oorawant  tharahy  laeura  no  raaponaiblllty^  nor  any 
ohllfatiOQ  ahataoavar)  and  tha  fact  that  tha  Oovam* 
■ant  nay  hava  fonulatad,  fomlahad,  or  in  any  way 
BuppUad  the  aald  dn«lnca>  apaeifleatlona>  or  othar 
data  la  not  to  ha  z«tardad  hy  laplleatlon  or  othar* 
vlaa  aa  In  any  Humar  lleanalnf  tha  holdar  or  any 
othar  paraon  or  eoxporatlon,  or  coovaylaf  any  rli^ta 
or  pamiaaloo  to  aaBufaetura,  uaa  or  aall  any 
patantad  Inrantion  tiiat  aay  In  any  way  ha  ralatad 
tharato. 


MATHEMATICS  RESEARCH  CENTER 


MATHEMATICS  RESEARCH  CENTER,  UNITED  STATES  ARMY 
THE  UNIVERSITY  OF  V^ISCONSIN 


Contract  No.:  DA-11-022-ORD-2059 


TRANSIENT  MAGNETOHYDRODYNAMIC 
FLOW  IN  AN  ANNULAR  CHANNEL 

i 

I  . 

I  M.  N.  L.  Narasimhan 

1 

I 

MRC  Technical  Summary  Report  #381 
February  1963 


Madison,  Wisconsin 


ABSTRACT 


Unsteady  parallel  flow  of  an  electrically  conducting  viscous 
Incompressible  fluid  In  an  annular  channel  In  the  presence  of  a 
radial  transverse  magnetic  field  Is  considered.  Assuming  the 
fluid  to  be  at  rest  at  the  Initial  moment,  the  velocity  distribution 
and  magnetic  field  components  are  obtained  In  terms  of  Bessel 
and  Lommel  functions  and  In  the  form  of  convolution  Integrals 
taking  the  longitudinal  pressure  as  an  arbitrary  function  of  time. 
Further,  taking  a  step  function  for  the  pressure  gradient,  these 
expressions  are  Integrated.  The  Influence  of  the  magnetic  field 
on  the  annular  flow  Is  Investigated.  Some  numerical  examples  are 


given. 


TRANSIENT  MAGNETOHYDRODYNAMIC  FLOW  IN  AN  ANNULAR  CHANNEL 


M.  N.  L.  Narasimhan 


Introduction. 

In  a  previous  paper  (Narasimhan  1963}  we  have  investigated  a  pulsating 
flow  in  an  annular  channel  with  a  radial  transverse  magnetic  field.  The  problem 
of  unsteady  magnetohydrodynamic  flow  in  a  rectangular  duct  has  been  solved 
by  Lundgren,  Atabek  and  Chang  (1961).  The  purpose  of  the  present  paper  is  to 
obtain  a  transient  solution  for  the  problem  of  unsteady  flow  in  an  annular 
channel  with  an  impressed  radial  transverse  magnetic  field.  We  consider  un 
infinitely  long  annular  channel  of  inner  radius  a  and  outer  radius  b,  shown 
in  Figure  1. 


FlgvAie  1.  Annular  channel  with  Impressed  radial  field. 
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We  consider  a  radial  magnetic  field  ^ ,  where  u  is  a  constant,  impressed 

across  the  channel.  In  practice,  it  is  possible  to  obtain  an  approximation  to 
the  desired  field  as  explained  by  Globe  (1959). 


Governing  equations. 

The  non-steady  flow  of  an  electrically  conducting  incompressible  fluid 
(mercury  for  Instance)  in  the  presence  of  a  magnetic  field  is  governed  by  the 
following  equations  (in  m.k.  s.  units) 


V  X  H  =  I  (1) 

VX  E  =  (2) 

V‘H  =  0  (3) 

■j  =  (r(E +(i  VX  H  )  (4) 

V* V  =  0  (^ 

p[^+(V-V)v]  = -Vp+  pv  V+pTxH  .  (6) 


The  first  four  equations  are  Maxwell's  equations,  which  govern  the  electromagnetic 
field,  Eq.  (5)  and  (6)  are  respectively  the  equation  of  continuity  of  fluid  and  the 
equation  of  motion.  In  the  above  equations  we  have  neglected  the  displacement 
current  and  assumed  that  the  permeability  and  dletectric  constant  are  the  same 
as  In  a  vacuum.  Also  we  have  neglected  the  free  charge.  We  use  cylindrical 
coordinates  (r,  6,  ?)  and  make  the  following  assumptions  in  the  manner  of 
Globe  (1959): 


a 

•^  =  0,  on  account  of  axial  symmetry. 

9  9 

v^.Ve  =  o  . 

Cl) 

We  assume  that  the  applied  field  ~  fixes  the  normal  component 

of  the  magnetic  field  at  r  =  a  and  r  =  b  for  all  values  of  z,  and 

that  this  is  the  only  field  Impressed.  Several  consequences  follow 

from  this  assumption. 

(a)  Now  from  (1)  and  (4)  we  have 

-  -  “6 

— t-  (..1) 

which  must  vanish  since  can  arise  only  from  an  applied  field 

E  or  free  charges  in  the  flow,  neither  of  which  exists,  and  because 

(VXHJ  must  vanish,  since  V  has  a  z-component  only.  Again 
z 

from  (1)  and  (4)  we  have 

-8H 

J^  =  alE^+fi(VXHy  =  -p<rV^HQ=  .  (a.  2) 

Thus  we  have 


Hence  from  (a.  2)  and  (a.  3)  we  get 

V  =-Ln5) 

z  (ji(r  f[z}  ’ 


(a.  4) 


a  function  of  z  alone. 
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Ftom  (5),  and  assumptions  (2)  and  (a.  4),  we  obtain 


=0 

dzUz)/  » 


1^2 

whose  solution  has  the  form  f(z)  =  Ae  ,  where  A  and  B  are 

constants.  This  would  make  V  a  constant  and  since  V  has 

z  z 

to  vanish  at  the  boundaries,  this  constant  would  be  zero.  Thus 

the  only  possibility  we  have  to  consider  is  that  f(z)  =  0,  so  that 

from  (a,  3),  H.  s  0  and  hence  j  also  vanish  everywhere  in  the  channel, 
o  r 

(b)  From  (3),  we  have 

dH  H  dH 

ar  ^  T  ^  dz  ■  °  » 

whose  solution  satisfying  the  boundary  conditions  can  be  obtained 
by  putting 

aH 

•7 

■  2F  ; 


az 


therefore 


aH 


ar 


H 

-  +  —  =  2F 


from  which  it  follows 


H  =  Fr  +  ^  , 

r  r  ^ 


r  r 

where  F  and  G  are  constants  at  any  Instant  of  time. 

Since  the  radial  magnetic  field  is  to  remain  at  its  prescribed 

(a)  (a> 

values  —  and  ~  at  r  =  a  and  r  =  b,  respectively,  for  all  values 

of  z,  we  obtain 

«  _  Q.  “  Q. 

-.Fa  +  ^. 

which  shows  that  F  =  0  and  G  =  w  . 
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Therefore 

8H 

a  0  and  *  0  or  rH^  =  constant. 

The  radial  field  must  then  be  equal  to  the  impressed  field  and  is 
unaffected  by  the  flow  and  is  independent  of  z  . 

Now  we  eliminate  J  and  E  from  (1),  (Z)  and  (4)  and  obtain 
the  following  equation 

»  -k  -k  1  7  -k 

(9) 

With  the  above  assumptions,  Equations  (6)  and  (9}  reduce  to  three  in  number: 

(10) 


^  -(i?-V)V  +  (V*V)H=j^V^H 


o  8H 


8v  8H 

_ z  uu _ z  8p  . 

P  at  "  r  8r  ■  ”  8z 


8H  8v  8^H  ,  8H 

_ Z  ucjq-  _ Z  _ Z  .  1  _ Z 

r  lr'~T+T~ 

Or 


(11) 


(12) 


and  v  are  functions  of  r  and  t  only.  It  follows  from  (11)  that 
z  z  8z 

must  be  independent  of  z  .  By  differentiating  (10)  with  respect  to  z,  it 
can  be  seen  that  ^  is  independent  of  r  also.  We  may  therefore  write 


8z 


P(t)  . 


Once  H  is  determined,  the  variation  of  p  across  the  channel  may  be 

Z 


found  by  integrating  (10): 

+  2  h!  =  “  PCt)  z  • 

Now  (11)  and  (12)  can  be  rewritten  as: 


8v  8H 

_ Z  UQ  _ Z 

'  8t 
8H 


8^ 


fir  ■  ‘’t'l  *  pv 


8r 


»  I  8v 

_Z  X  2 
2  r  8r 


,  8v  8^H 

2  _  UMO-  Z 


-z  1  ®“z 

8t  "  r  8r  "  2  r  8r 

Ot 


(13) 


(14) 
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The  initial  conditions  are: 


V  (r,t)  =  0  ,  at  t  =  0  .  (15) 

z 

H  Cr,t)  =  0  ,  at  t  =  0  .  (16) 

z 

The  boundary  conditions  of  the  problem  are: 


v^(r,t)  =  0 

at  r  =  a 

(17) 

v^(r,t)  =  0 

at  r  3  b 

(18) 

H^jCr,t)  =  0 

at  r  =  b 

(19) 

dH 

— = « 

at  r  =  b  . 

(20) 

(17)  and  (18)  are  the  no  slip  conditions. 

(19)  follows  from  the  fact  that  J  has  a  e-component  only  so  that  the  currents 
in  the  annular  channel  are  like  those  in  an  infinite  solenoid.  These  currents 
produce  no  field  for  r  >  b  and  since  there  is  no  impressed  field  in  the 
z-dlrectlon,  continuity  of  the  tangential  component  of  H  requires  (19)  to  be 
true.  (20)  is  Justified  as  follows: 


but 


"  "  "air 


jg  =  (r|*(VXH)g 

.  Hence  J.  must  also  vanish  there  and 

P 


dH 
_ z 

8r 


and  V  must  vanish  at  r  =  b 


too. 


#381 


7 


Now  we  non-dimenelonallze  the  equations  (13)  to  (20)  by  introducing  the 
following  dimensionless  quantities: 


Let 


(v/a)  • 


tv 

a  (hktv)* 
H 


H 


(v^){(rpv) 


♦(t)=  -a^p"^ 


2  2 


p2  ^  ^  Y  »  and  let  I  =  6  ; 


p  is  a  farm  of  the  Hartmann  number. 


When  these  non-<limensional  quantities  are  introduced  into  the  equations  (13) 
to  (20),  the  latter  take  the  form: 


LSI  ,^.LS1.S.SS. 

Y  8t  2  V  8X  \  8\ 


♦(t) 


(21) 


8X 


^  IB.Oi  .XM  .J.8JL 

^  8t  "  2  X  8X  X  8X 

9\ 


(22) 


with  the  conditions: 


V(X,t)  =  0,  at  t  =  0  (23) 

H(X,t)  =  0,  at  t  =  0  (24) 

V{1,T)  =  0,  (25) 

V(6,t)  -  0  ,  (26) 
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•  8** 


H(8,t)-0  ,  (27) 

1^  (K,t)  »  0  at  \  «  6  .  (28) 


TM  po^Utton. 


First  let  us  eliminate  t  from  the  equations  by  taking  a  Laplace  transform 
with  respect  to  t  .  Let  s  be  the  transformed  parameter  and  V ,  H  and  ^ 
denote  transformed  variables.  Thus  our  problem  reduces  to: 


8^V  ^1  W  6  8H 

^.2*  \  8K  K  8K  “ 
oK 


(29y 


afs  1 


(30) 


with  the  conditions: 

V(l,s)  =  0  , 
V(6,8)  =  0  , 
H(6,8)  =  0  , 


8H,.  , 

8\ 


(31) 

(32) 

(33) 

at  X  s  6  .  (34) 


The  solution  of  this  system  of  equations  for  general  values  of  y  s  (ixo-  v)^ 
runs  into  difficulties.  Moreover,  in  actual  physical  situations,  it  will  be 

sufficient  to  solve  the  system  for  y«  1,  since  for  most  of  the  incompressible, 
electrically  conducting  viscous  fluids  on  the  surface  of  tl.i  earth,  y  <<  1  • 


#381 
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For  instance  (Lundgren,  Atabek,  Chang  1961): 


Table  1 


Fluid 

V  =  (u<rv)2 

Hg  (20'C) 

3.56  Xio"* 

Na  (500“C) 

-3 

1.55  XIO 

-3 

Pb  (500*0) 

0.47  X  10 

Hence  we  shall  assume  y  <<  1  and  s  •  The  equations  (29)  and  (30)  , 
after  neglecting  terms  of  order  ,  reduce  to: 


(35) 


and 


afj_  1  ^ 


X  av 


«  0 


(36) 


From  (36)  after  integrating  once  and  using  (32)  and  (^34)  we  obtain 


an 

ax 


(37) 


Substitution  of  (37)  into  (35)  yields  the  following  non-homogeneous  modified 
Bessel  equation: 


2  a^v 


ax 


+  x^- 

2  ^  ^  ax 


(x^^  +  p2)v  =  -x^; 


(38) 
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The  general  solution  of  (38)  is  therefore  of  the  form 


V(X  ,  s)  =  Aj  Ip(X  *s/177 )  +  A2Kp(XN/  s/y  )  +  J  s)  p(iV  ^  s/y  )  ,  (39) 


where  the  functions  1^,  are  modified  Bessel  functions  of  the  first  and  second 
kind  respectively  and  of  order  f3,  and  s^  is  a  Lommel  function  (Watson, 

1944)  which  occurs  as  a  particular  Integral  of  (38)  .  A  and  A  are  arbitrary 

i  M 

constants.  This  particular  integral  can  be  written  either  as  ascending  or 
descending  power  series  in  the  argument  (Watson  1944)  accc^rding  as 
l*p#-(2p-l)  or  l*p=-(2p-l),  where  p  is  a  positive  Integer.  For  the 
sake  of  simplicity  we  shall  demonstrate  the  solution  when  1  *  -(2p-  1)  . 

This  condition  is  satisfied  when  p  =  an  odd  positive  integer  or  zero.  Thus  we  shall 
demonstrate  the  solution  when  p,  the  Hartmann  number  is  an  odd  positive 
Integer  or  zero. 

Hence  in  this  case  we  obtain  the  complete  solution  of  (38)  as: 


V(\  ,s)« 


ix 

s 


Sj^p(i6Vs/Y  )  {Kp(N/s/Y)Ip(\N/s/Y)  -Ip(^/s/Y)Kp(X  ^  s/y)} 

-S^  p(i^/^){Kp(6^/s/Y  )  Ip(XN/s/Y  )-I^6Vs/y 

)Kp(xVs/Y)} 

+  Sj^p(i\N/^) 

Ip('^)Kp(6^/I77  )  -  )  Kp('^^  ) 

(40) 


#381 


-11- 


where 


(ml  +  1)^-v^ 


{ttx+l)^-v'^>{(>i+3)^-v^ 


+  •  .  • 


(41) 


and  (i  *  V#  -(2p-l),  p  being  a  positive  integer. 

Now  taking  the  inverse  Laplace  transform  of  (40)  and  using  the  convolution 
theorem  (Carslaw  and  Jaeger  1941)  we  obtain 


T 

V(V,t)  =  2  a  /«|►(T-|)exp(-Ytt„l)d4  , 

n=l  "  0  " 

where 


(42) 


A 

n 


(V“n‘>  VV>- W>  W» 


(4  3) 


Jp  and  Yp  being  Bessel  functions  of  the  first  and  kind  respectively  and  a  's 
(n  =  1, 2, . , .  00 )  are  the  roots  of  the  equation 

Jp(4  Yp(«S)  -  Jp(a6)  Yp(a)  =  0  ,  (44) 

P  being  a  positive  odd  integer.  These  roots  a  's,  n  =  l,2,.  ..oo  are  all  known 
to  be  real  and  simple  since  6  and  p  are  real.  [Carslaw,  Conduction  of  heat 
(London  1922)  p.  128]. 
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V(V,t) 


HCV,t) 


To  Integrate  (42)  we  have  to  specify  time  dependence  of  the  pressure  gradient. 

Step  funcUon  pressure  gradient. 

Let  us  take  ^r)  as  a  step  function.  That  is  ,  let 


♦(T)  = 


0  for  T<  0 


for  T  >  0 


(45) 


Then  after  integration.  Equation  (42}  gives  us 


1  "Y«  T 

1-e  '  n 

I 

(46) 


By  taking  the  inverse  transform  of  (37)  and  using  (26)  and  (28)  we  have 


m  ■  _  tT) 
8V  ■  X 


i47) 


If  we  now  substitute  (46)  into  (47)  and  integrate,  we  obtain  after  using  (27): 


“n  2  2  2 

n=l 


,  “va  T 

1-e  ^  n 


(48) 
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where 

-  u^(6)  -  -  p  /  d\ 

6 

VY  (Xa  ) 

w^CX>)-Wn(6)  .-p/  P I  dX  . 

6 

Thus  (46]  and  (48)  give  the  velocity  field  and  magnetic  field  respectively. 

Discussion  of  the  solution. 

The  expressions  (46)  and  (48)  for  V  and  H  contain  a  Ume-lndependent 
part  and  a  time-dependent  part.  The  first  parts,  namely  the  time-independent  parts 
represent  the  fully  developed  values  of  V  and  H,  respectively. 

In  the  time-dependent  parts  which  contain  exponential  terms,  we  find  that 
starting  from  rest  V  and  H  grow  to  their  developed  values  without  ever 
exceeding  them.  Since  ci^*s  are  all  real,  it  is  found  that  periodic  or  even 
partially  periodic  flows  are  Impossible  with  the  above  choice  of  step  function 
pressure  gradient  unlike  the  case  of  a  similar  type  of  flow  in  a  rectangular  duct 
with  a  transverse  magnetic  field  (Lundgren,  Atabek  and  Chang  1961).  Thus  it  is 
found  that  with  a  suddenly  applied  pressure  gradient  the  flow  of  conducting 
incompressible  fluids  under  a  radial  transverse  magnetic  field  is  damped  and 
partially  periodic  flows  are  impossible.  We  have  calculated  numerically  the 
velocity  distribution  when  the  radius  ratio  of  the  annulus  -  =  2  .  The 

d 

hydrodynamic  and  hydromagnetlc  cases  corresponding  to  the  values  of  the  Hartmann 
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number  p  s  o  and  p  a  5  have  been  obtained  for  purposes  of  comparison.  The 
roots  or^'s  of  the  equation 


Jp(«)Yp{«4-Jp{ar6)Yp(or)  =  0;  6-^-. 


=  r  =  2 


have  been  calculated  by  using  the  formula  for  a^'s  (Gray  and  Mathews  1895): 

"  "  k* 


where 


K  -  OIL  a  -  4fm-lWm-25U6^-ll 


3(86)-' (6-1) 


r  =  32Lm-l)(m^-114m4l073){6^-n  ^  ,  4p2  . 

5(86)®  (6-1) 


the  first  ten  roots  are  tabulated  for  p  s  0  and  p  =  5  in  table  2  . 

The  Lommel  function  s^  ^x)  for  p  ■  0  reduces  to  the  well  known  relation 
(Watson  1944), 


where  JJx)  is  the  Bessel  function  of  the  first  kind  and  order  zero  and  s  (x) 

u  1,0'  ^ 

has  been  computed  thus.  The  Lommel  function  ^(x)  f or  p  =  5  has  been 
computed  from  equation  (41)  and  for  large  arguments  x  the  asymptotio  relation 
(Watson  1944) 


s  (x)  s  x' 


has  been  used. 


#381 


15- 


The  results  of  the  computation  are  given  by  the  velocity  profiles  in 
Figure  2  and  Figure  i,  obtained  as  functions  of  the  radial  coordinate  K  .  In 
Figure  2  the  velocity  profiles  for  the  hydrodynamic  flow  «  0),  have  been 
calculated  when  a  o.l,  0.25,  0. 5,  1.0  and  w  .  In  Figure  3,  the  velocity 
profiles  for  the  hydromagnetio  flow  XP  °  3),  have  been  calculated  when 
yr  3  0. 1,  0. 5  and  «  ,  From  these  figures,  it  is  found  that  the  effect  of  the 
magnetic  field  is  to  flatten  the  velocity  distribution  and  to  shorten  the 
development  time.  In  Figure  3,  it  is  clearly  seen  that  the  flow  development 
time  for  the  hydromagnetlc  flow  is  very  much  less  compared  to  the  development 
time  for  the  hydrodynamic  flow. 

Summary 

In  this  paper  the  unsteady  flow  of  a  conducting  fluid  through  an  annular 
channel  under  a  radial  transverse  magnetic  field  and  a  suddenly  applied  pressure 
gradient  has  been  considered.  For  most  conducting  incompressible  fluids  for 
which  Y  ■  (tA<rv)^<<  1,  the  response  of  the  fluid  is  damped  and  partially  periodic 
flows  are  found  to  be  impossible.  The  effect  of  the  magnetic  field  on  the  fluid 
flow  is  found  to  flatten  the  velocity  distribution  and  to  shorten  the  flow 
development  time. 
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The  first  ten  roots  ar^'s  of  the  Equations  Jp(a)  Yp(ar  6)  -  Jp(«6)  Yp(«)  »  0  for 


I  thank  Professor  D.  Greenspan  for  making  available  some  assistance  in 
numerical  caloulations.  My  thanks  are  also  due  to  Mr.  Donald  Van  Egeren  for 
performing  some  of  the  numerical  calculations. 
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